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Abstract 



The supersymmetric flavor, CP and Polonyi problems are hints that the 
fundamental scale of the soft supersymmetry breaking parameters may be 
above a TeV, in apparent conflict with naturalness. We consider the possibil- 
ity that multi-TeV scalar masses are generated by Planck- or unification-scale 
physics, and find the conditions under which the masses of scalars with large 
Yukawa couplings are driven, radiatively and asymptotically, to the weak scale 
through renormalization group evolution. Light third generation scalars then 
satisfy naturalness, while first and second generation scalars remain heavy to 
satisfy experimental constraints. We find that this mechanism is beautifully 
realized in the context of grand unified theories. In particular, the existence 
of right-handed neutrinos plays an important role in allowing remarkably sim- 
ple scenarios. For example, for SO(10) boundary conditions with the squared 
masses of Higgs scalars double those of sleptons and squarks, we find that the 
entire scalar mass scale may be increased to 4 TeV at the unification scale 
without sacrificing naturalness. 
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I. INTRODUCTION 



The Achilles heel of supersymmetry lies, arguably, in the difficulty of satisfying natural- 
ness while simultaneously decoupling supersymmetric effects from low-energy experimental 
observables. Naturalness is typically taken to require supersymmetric particle masses below 



y (4:7i/a)Mw ~ 1 TeV, a few times the weak scale. In contrast, without an understanding 
of small scalar mass splittings, mixing angles, and CP- violating phases, a variety of ex- 
perimental measurements suggest that supersymmetric scalar masses must lie well above a 



The experimental constraints have varying degrees of significance and model dependence. 
Examples include the following: 

• The kaon system. Constraints from the K system depend on the flavor and chirality 
structure of the scalar quarks. Roughly speaking, however, for moderate degeneracies, 
the constraint from €k is most severe; it requires squark masses of order 100 TeV. 
Likewise, the K — K mass splitting needs squark masses of order 10 TeV. (See Ref. |I[ 
for details.) 

• Electric dipole moments. The electric dipole moments of the neutron and electron 
typically require multi-TeV squark and slepton masses if the relevant phases are not 
suppressed [[|. EDMs are flavor-conserving and so cannot be suppressed by scalar 
degeneracy; for example, they imply severe constraints on scalar masses (or phases) 
even in models with gauge-mediated supersymmetry breaking. ||. 

• The proton lifetime. Lower bounds on the proton lifetime place strong constraints on 
supersymmetric grand unified theories (GUTs), especially for large tan/3 The pre- 
dicted decay rates depend on the (super)heavy and the light supersymmetric particle 
spectra; they are greatly suppressed for multi-TeV scalar masses. 

• The Polonyi problem. Many supersymmetric theories contain late decaying scalar fields 
that ruin the successful predictions of Big Bang nucleosynthesis 0. This problem 
is independent of scalar degeneracies or phases, but may be solved if the scale of 
supersymmetry breaking is of order 10 TeV 0. 

In addition, many models predict scalar masses hierarchically larger than gaugino 
masses J7J. Given the current bounds on gaugino masses, scalar masses must then be well 
above a TeV. 

The conflict between naturalness and experimental constraints may be resolved by ob- 
serving that, roughly speaking, naturalness restricts the masses of scalars with large Yukawa 
couplings, while experiment constrains the masses of scalars with small Yukawas ||. Nat- 
uralness affects particles that are strongly coupled to the Higgs sector, while experimental 
constraints are strongest in sectors with light fermions that are plentifully produced. This 
suggests that naturalness and experimental constraints may be simultaneously satisfied by 
an "inverted hierarchy" approach, in which light fermions have heavy superpartners, and 
vice versa. In particular, third generation scalars with masses ~ 1 TeV satisfy nat- 

uralness constraints, while first and second generation scalars at some much higher scale 




TeV. 
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"^heavy avoid many experimental difficulties. A number of possibilities have been proposed 
to dynamically generate scalar masses at two hierarchically separated scales |p|-|TT|]. 

In this paper, we shall investigate a new mechanism for generating an inverted hierarchy. 
We will assume that all scalar masses are of order m heavy when produced at some scale, say 
the GUT or Planck scale. We will then use renormalization group evolution to create an 
inverted hierarchy at the weak scale. This mechanism works because the renormalization 
group equations automatically suppress the masses of scalars associated with fermions which 
have large Yukawa couplings. The masses of these scalars (and only these scalars) are driven 
to miight- This mechanism was investigated recently in models with large top and bottom 
Yukawa couplings (and all others small) |12|]. It was shown that, for particular high scale 
boundary conditions, all scalars with large Yukawas may be driven asymptotically to mu g ht 
in the infrared. Stops, sbottoms, and Higgs scalars are then naturally light at the weak 
scale, while the masses of all other scalars remain heavy. In this way a very large hierarchy 
can be created dynamically. 

In what follows we will consider the well-motivated supersymmetric GUTs, with special 
emphasis on SO(10) unified theories. As is well-known, such theories, among other merits, 
successfully predict gauge coupling unification and naturally accommodate massive neutri- 
nos. In the present context, however, they have even more virtues. First, they naturally 
provide the universal and large third generation Yukawa couplings required to implement 
the inverted hierarchy. Second, as we will see, the existence of a right-handed neutrino 
plays an important part in creating a plausible fixed point structure. These virtues make 
grand unified theories the natural setting for the radiatively generated inverted hierarchy 
mechanism. 

We will begin in Sec. |I| with a discussion of the fixed point structure of the scalar 
mass renormalization group equations (RGEs). In Sec. [ITI| we will analyze the inverted 
hierarchy in grand unified models with right-handed neutrinos. Finally, in Sec. [IV] we will 
consider inverted hierarchies generated in supersymmetric GUTs with evolution between the 
Planck and GUT scales. In each case we illustrate the size of the hierarchies that may be 
achieved using the scalar mass fixed-point framework. We shall see that the scalar masses 
may be pushed as high as ~ 4 TeV. Such a large hierarchy significantly ameliorates most 
experimental difficulties. More detailed conclusions, as well as a discussion of implications 
for low energy experiments and high energy colliders, are contained in Sec. [V]. 

In the Appendix we collect the RGEs used in this analysis. For models above the GUT 
scale, these RGEs include, for the first time, the two-loop pure Yukawa terms. They also 
correct some one-loop results previously given in the literature. 



II. SCALAR MASS FIXED POINTS 

Supersymmetric theories are natural if the conditions for electroweak symmetry breaking, 

1 2 _ m 2 Hd - m? Hu tan 2 g 2 
2 mz ~ tan 2 /3-l // ' 

2m 2 = (m 2 + m 2 + 2fi 2 ) sin 2(3 , (1) 
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are free of large cancellations. [] These conditions apply to the RGE-improved effective po- 
tential, so the supersymmetric parameters of Eq. (|l|) must be evaluated at the weak scale. 
A priori, naturalness bounds only the parameters entering Eq. ([!]), and, for example, even 
the requirement on m# d is relaxed for large tan f3. Nevertheless, naturalness affects all other 
supersymmetric parameters to the extent that they renormalize these relations. 

For all the reasons given in the previous section, we are motivated to consider theories 
in which the scalar masses are of order m^y at some high scale Mo, which we take to be 
the GUT or Planck scale. For arbitrary initial conditions, it is clear that the scalar masses, 
including mn u and are of order mh ea vy at all scales, and the theory must be fine-tuned. 
However, if the RGEs have (approximate) infrared fixed points where these scalar masses 
vanish, appropriate boundary conditions will lead to weak-scale masses that are small and 
insensitive to mheavy, thereby preserving naturalness. 

We must therefore identify the infrared fixed points of the scalar mass renormalization 
group equations. The RGEs for supersymmetric theories are well-known (up to two- loops), 
but complicated However, in the scenarios we are considering, we can make several 



simplifications and easily extract the essential features of the RGEs. 

In a general supersymmetric theory, the one-loop RGEs for scalar squared masses m 2 
are schematically of the form 

m 2 Ym 2 + aM 2 /2 - Y A 2 , (2) 

where M1/2 and A represent gaugino masses and trilinear scalar couplings, respectively.^ 
Here, and in the rest of the paper, we use the following notation and conventions: t = 
In (Mq/Q 2 ), ( ) = j- t , a = g 2 /16ir 2 = a/4ir, and Y = h 2 /lQiT 2 , where g and h are gauge and 
Yukawa couplings, respectively. For each term in Eq. (0), overall signs are determined as 
shown, but (positive) numerical coefficients are suppressed. 

From Eq. (fj) we see that large M\/2 and A parameters typically destabilize light scalar 
masses. Therefore we require^ M\/2,A ~ mi ig ht at the scale M . Note that the hierarchy 
Mi/2, A<m at M is natural and may be generated in a number of ways, for example, by 
an approximate .R-symmetry, or through mechanisms that differentiate dimension-one and 



dimension-two soft supersymmetry breaking parameters. Examples are given in Refs. [jl2| , pT4 

Let us also assume that all large Yukawa couplings are unified at Mq. Then, if we ignore 
differences in the Yukawa evolution, we find that the RGEs for the scalar masses take the 
simple form 



x The parameters mn u and mn i are the soft supersymmetry breaking Higgs boson masses, m 2 is 
the soft bilinear scalar coupling of the two Higgs doublets, /i is the Higgsino mass parameter, and 
tan/3 = (H®)/(Hd) is the usual ratio of Higgs vacuum expectation values. 

2 We have omitted a hypercharge trace term which vanishes for supersymmetric GUTs. 

3 Gaugino masses satisfy (^Mi/ 2 /ckj = at one-loop; they cannot be driven from mheavy at Mq 
to TTihght in the infrared. For ^-parameters, the story is more complicated. We will simply assume 
A ~ Mi j2 ~ flight at all scales. 
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mj 



-YN ijm j 



(3) 



where the subscripts run over all scalar fields in the theory, and N is a matrix of positive 
constants determined by color and SU(2) factors. This set of equations is easily solved by 
decomposing arbitrary initial conditions into components parallel to the eigenvectors of N, 
each of which then evolve independently. Indeed, if N has eigenvectors e$ with eigenvalues 



Aj, the initial conditions 



m 2 (t = 0) = J2 C ^ 



evolve to 



m 2 (t = t f ) =J2 °i e i ex P ~\ / 

L <J 



rtf 

-Xi I Ydt 



(4) 



(5) 



We see that eigenvectors with large positive eigenvalues are asymptotically and rapidly 
crunched to zero. For initial conditions dominated by such eigenvectors, the final masses 
are greatly suppressed relative to their initial values. 

Of course, the size of the scalar mass hierarchy depends not only on the initial conditions, 
but also on the evolution interval tf and the initial value of the Yukawa coupling. We expect 
a large hierarchy for large tf and Yukawa couplings near their quasi- fixed points [15|]. In the 



following sections, we determine numerically the suppression factors, or crunch factors, that 
are possible in supersymmetric GUTs both below and above the unification scale. 



III. MSSM WITH RIGHT-HANDED NEUTRINO 

In this section, we consider a model with the particle content of the minimal super- 
symmetric standard model (MSSM) with a right-handed neutrino N. We assume that soft 
supersymmetry breaking masses are generated at the unification scale, as may be the case, 
for example, in strongly-coupled string theories where the GUT scale is also the string and, 
hence, the supergravity scale. 

The superpotential is given by 

W = h t H u QU + h b H d QD + h T H d LE + h n H u LN , (6) 

where the matter fields Q, U, D, L, E, and N are those of the third generation and all other 
Yukawa couplings may be neglected for our analysis. We also neglect off-diagonal scalar 
masses for the moment; their effects and constraints will be discussed in Sec. |V|. 

For this model, omitting gaugino masses and A terms as discussed in Sec. |J, the RGEs 
for scalar masses are 
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where m 2 = ijn 2 Hu) mfj, uiq, m 2 D , m 2 Hd , m 2 L) m 2 El m 2 N ) T . At the GUT scale Mq, the Yukawa 
couplings are unified with Y = Y t = Yj, = Y T = Y n . If we assume, for the moment, that they 
remain approximately degenerate at lower scales, the evolution equations simplify to 



m 



-YNm 2 



(8) 



where 



N 



4 3 3 

2 2 2 
112 

2 
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10 





10 1 
10 1 
10 
2 
4 110 
12 11 
2 2 2 
20000202 



(9) 



The eigenvectors of N, and their associated eigenvalues Aj, are 

ei: 8, (2, 1,1, 1,2, 1,1,1) 
e 2 : 6, (2, 1,0, -1,-2, 0,-1,1) 
e 3 : 4, (0,1, 1,1, 0,-3, -3, -3) 
e 4 : 2,(0,-1,0,1,0,0,-3,3) , 



(10) 



along with four eigenvectors with zero eigenvalue. Components of the initial conditions along 
e l5 and to a lesser extent e 2 , e 3 , and e 4 , are rapidly suppressed during renormalization group 
evolution by factors of exp[— Aj / Ydt]. 

Note that the boundary condition given by ex, 



rrir 



m, 



m 



D 



m, 



m% 



m 



N 



-m 



(11) 



is remarkably simple, and is consistent with minimal SO(10) unification, in which all matter 
fields arise from a single 16 multiplet, both Higgs fields are contained in a single 10 mul- 
tiplet, and the GUT breaking SO(10) -> SU(3) x SU(2) x U(l) takes place in one stepfj 
This simplicity is not automatic. For example, equivalent analyses without a right-handed 
neutrino yield eigenvectors that are far from simple, and moreover, are not compatible with 
GUT unification. The right-handed neutrino plays a critical role in leading us to a plausibly 
simple boundary condition. 

In reaching the above conclusions, we have made a number of approximations that must 
be examined. First, we have neglected the fact that the Yukawa couplings evolve indepen- 
dently, and, in particular, the fact that the leptonic and hadronic Yukawas differ significantly 



4 Of course, such a minimal unification scenario is not consistent with the observed light fermion 
spectrum and CKM mixing. However, these issues may be resolved, for example, by non- 
renormalizable couplings, which are irrelevant to our analysis of the large 0(1) Yukawa couplings. 
The minimal model may therefore be used to illustrate our results. 
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at lower energy scales. Second, we have not taken into account the decoupling of the right- 
handed neutrino iV at some scale. Both of these effects imply that the above eigenvectors 
do not remain eigenvectors during the full renormalization group evolution. Finally, as we 
will see below, there may also be corrections from numerically significant two-loop terms. 

To determine the importance of these effects, we now solve the RGEs numerically. It is, 
of course, possible to include all terms up to two-loops in the numerical analysis. However, 
doing so requires specifying many additional parameters, and the analysis becomes highly 
model-dependent. Fortunately, many terms have only a small effect on the overall fixed point 
behavior and therefore may be neglected. In what follows, we will systematically omit all 
terms of order m 2 ight in the m 2 RGEs; this removes all dependence on the unknown gaugino 
mass and A parameters. For consistency, we also omit all other terms of similar size, and 
keep only the one-loop and leading two-loop terms. The resulting RGEs are schematically 
of the form 

a = ba 2 + a 2 Y 

Y = -Y 2 + aY + Y 3 (12) 
m 2 = —Ym 2 + Y 2 m 2 , 

where h is the one-loop /3-function coefficient. The complete RGEs up to this order are 
presented in the Appendix. 

The approximation used is formally equivalent to keeping the leading and next-to-leading 
order terms in an expansion in 0(V&,Y). This expansion is motivated by the fact that, at 
the unification scale, Yukawa couplings may be very large near their quasi-fixed points. For 
example, for h ~ 2, we have Y ~ 1/40, while a ~ 1/300. The dominant two-loop terms 
are therefore those of Eq. (0). We will see that these terms may be significant. On the 
other hand, while possibly significant, they should not be too large: large two- loop effects 
signal a breakdown of perturbativity, and our analysis is unreliable in these regions. We will 
comment on such parameter regions below. 

With these approximations, the theory is completely specified by two parameters, ha, 
the universal third generation Yukawa coupling at the GUT scale, and m^r, the right-handed 
neutrino mass. In Fig. [I], we show the renormalization group evolution of the scalar mass 
parameters from the GUT scale Mq — 2 x 10 16 GeV to the weak scale, for ha = 2 and 
m^r = 10 13 GeV. We find that large scalar mass suppressions are possible. Indeed, even 
^heavy ~ 4 TeV can give rise to third generation scalars with masses m light < 1 TeV. Note 
also that the hierarchy is created rapidly, in the first few decades of evolution. This is the 
region in which the gauge interactions, and hence, the Yukawa couplings, remain roughly 
universal. Below this scale, the Yukawa couplings split substantially, but by then, the 
hierarchy is already created and cannot be destroyed. The assumption of universal Yukawa 
couplings throughout the evolution interval is therefore a convenient fiction. 

In our analysis, we have neglected effects of order mii g ht, which clearly have little impact 
on the overall suppression of the scalar mass scale. They do, however, modify the RGEs 



5 Note that a 2 ?ri^ eavy terms have been omitted in the m 2 RGE. However, given the hierarchies we 
are able to achieve (see below), these terms are still only of order m 2 ight . 
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FIG. 1. The renormalization group evolution of the Higgs (dashed) and third generation squark 
(solid) and slepton (dotted) squared masses in the MSSM with a right-handed neutrino for the 
boundary conditions of Eq. (|l~l] ) with \iq = 2 and rajy = 10 13 GeV. First and second genera- 
tion scalar masses (not shown) are approximately renormalization group invariant. At the weak 
scale, m L = 1270 GeV, m E = 830 GeV, m Hu = 710 GeV, m D = 690 GeV, m Q = 570 GeV, 
mu = 420 GeV, and rriH d = 50 GeV. Note, however, that neglected effects of order rnf ight modify 
solutions in the shaded region. The suppression factor for this case is S = 20 (see text). 



at the weak scale when the scalar masses are of order m^t- By omitting such effects, we 
forfeit the possibility of investigating a number of topics, including the details of electroweak 
symmetry breaking, as well as quantitative determinations of finite corrections to fermion 
masses and the residual fine-tuning required to keep the Z mass light. Here, we note only 
that our scenario shares the problem of proper electroweak symmetry breaking generic to 
all supersymmetric scenarios with large tan (3 | I6| . Several possible solutions have been 
discussed in the literature. For example, if the GUT group is broken not in one step, but 
in a chain beginning with SO(10) — > SU(5) x U(l), the resulting U(l) D-terms may help 
induce proper electroweak symmetry breaking |l7j . These terms are parametricaliy of order 



m Lavy) but are suppressed by the large U(l)x charge of the singlet. More generally, another 
source for C(mf ight ) (but not C'(^heavy)) perturbations could be assumed. 

The large radiatively-generated hierarchy evident in Fig. [I] exists for a wide range of 
parameters and ho- To quantify the hierarchy, we characterize the scale of a given 
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FIG. 2. The suppression factor S denned in the text for the MSSM with a right-handed neutrino 



N with initial boundary conditions given by Eq. (11). The parameter is the scale at which 
the right-handed neutrino decouples, and Kq is the value of the universal Yukawa coupling at the 
GUT scale M G ~ 2 x 10 16 . 



scalar spectrum at renormalization scale Q by the quantity 

m 2 (Q) =Av[\m 2 (Q)\] , (13) 

where the average is taken over all scalar fields in the theory, properly weighted by color and 
SU(2) factors. (This definition of rh 2 (Q) is invariant under possible shifts from D-terms.) 
We then define the suppression, or crunching, factor 

8m^r, (14) 
m z {mw) 

where the weak scale m\y is taken to be 1 TeV. In Fig. |] we plot S in the (mjv, ha) plane. 
For a broad range of h t near its quasi- fixed point, we see that S ~ 20 is possible. For 
smaller values of h t , the scalar masses do not approach their fixed point as quickly, and for 
larger values, two- loop effects become important. Indeed, the two-loop terms give positive 
contributions to Y, and therefore increase the Yukawa suppression of scalar masses. However, 
they also give positive contributions directly to these masses. (See Eq. fll2]).) The latter 
effects dominate, so the two-loop terms degrade the hierarchy. For ht larger than shown in 
Fig. 0, the two-loop effects become so large that perturbation theory cannot be trusted. 

The definition of S is, of course, somewhat arbitrary, but its size gives an indication of 
the reductions in the scalar masses that are possible from radiative effects. In particular, 
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S provides a rough measure of the reduction of the scalar contribution to fine-tuning and 
is a measure of mheavy/^iight- From Fig. || we see that, if values of of order 1 

TeV are considered natural, values of m heavy up to ~ 4 TeV are acceptable. This is then the 
scale of the matter fields of the first and second generation, which, given their small Yukawa 
couplings, do not renormalize significantly. 

A scale of 4 TeV does not completely solve the most stringent flavor and CP problems 
related to the K system, but it does reduce them significantly with respect to the typical 
case with squark masses of 1 TeV. In addition, several of the other experimental difficulties 
mentioned in Sec. | are also resolved. Detecting scalars with such heavy masses will pose a 
serious challenge to future colliders. 

Finally, we note that the values of S are largely independent of mjv- For < 10 13 GeV, 
the hierarchy is already created by the time the neutrino decouples, and so is insensitive to 
the decoupling scale. However, even for approaching Mq, a large hierarchy is possible. 
This may be understood as follows. When N decouples, the new evolution matrix N' is the 
upper 7x7 block of AT" with Y n = 0. The matrix N' has 3 positive and 4 zero eigenvalues. 
In terms of cos 3$ = 17/64, the 2 largest eigenvalues are 3 (cos 9±y/3sm9 + 16) w 7.5,5.6; 
their corresponding eigenvectors e' x and e' 2 are also quickly damped. If we decompose the 7- 
dimensional truncation of e± into these new eigenvectors, ef = (2, 1, 1, 1, 2, 1, 1) = Sj =1 c[e[, 
we find that the decomposition is dominated by e[ and e' 2 components: for the remaining 
components, | J2j=3 c i e i\/\ et i\ = 0.20. Thus, the decoupling of N perturbs the fixed point 
structure only slightly — when iV decouples, the scalar mass spectrum projects mainly on to 
the new eigenvectors with largest eigenvalues, and the rapid exponential damping continues. 
Note that the special case m N = M G is equivalent to the case of an SU(5) GUT with 
universal Yukawa couplings and minimal SU(5) particle content. The exact eigenvectors for 
the SU(5) fixed point system are complicated. However, the above analysis shows that e^ r 
is nearly a linear combination of eigenvectors with large eigenvalues, so large hierarchies can 
be generated from simple boundary conditions in the SU(5) well. 

IV. GUTS ABOVE THE UNIFICATION SCALE 

An inverse hierarchy may also be generated in supergravity theories through renormal- 
ization group evolution above the unification scale. In this section we consider models in 
which scalar masses are generated at the supergravity scale, which we take to be the reduced 
Planck scale M* = 2.4 x 10 18 GeV, and then evolved to the GUT scale M G ~ 2 x 10 16 GeV. 
Even though the evolution interval is only two decades, we again find that substantial hier- 
archies may be created for simple initial scalar spectra. We will consider two generic cases 
for illustration. 

A. SO(10) 

We begin by considering SO(10) models, that is, the class of models with superpotential 

w = h t ^n 2 + . . . , (15) 

where the matter fields are contained in ip (16) and both up-type and down-type Higgs 
fields are contained in a single multiplet: 7i 2 (10) ~D iH u ,Hd), where H u , Hd are 5 and 
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5 representations of the SU(5) subgroup, respectively. We allow arbitrary additional field 
content, subject only to the constraint that additional superpotential couplings are small 
relative to h t . Note that two Higgs models, with superpotential W = h t ipipT-C2 + h b ipipT-ti + . . . 
are, for the purposes of our analysis, equivalent to a one Higgs model in the limit h t = h b 
(see Appendix). To the extent that the couplings are unified in this more general class of 
models, this analysis also applies. 

The RGEs for this class of models are given in the Appendix. Note that additional 
interactions and an extended Higgs sector are required to break SO(10) and generate fla- 
vor structure. In general, this leads to a set of RGEs that are highly model-dependent. 
Fortunately, though, in the limit we consider, the RGEs simplify considerably. First, we 
assume that all additional Yukawa couplings are small, so their impact on the RGEs may be 
neglected. Second, the presence of additional fields with standard model quantum numbers 
modifies all RGE terms with two or more powers of a. Fortunately, as argued in Sec. |TJ, 
the two-loop terms of this form are highly suppressed. (We have checked that they are 
insignificant for reasonable field content.) Third, extra fields with standard model quantum 
numbers also change the one- loop /3-function coefficient 610. To good approximation, this 
is the only important effect, so the RGEs of these models may be studied simply by taking 
610 as an arbitrary free parameter. We will consider the range —3 < 610 < 20, for which 
1/26 < a(M0 < 1/8. 



We now proceed as in Sec. |HI| . The one-loop scalar mass evolution is given by rh 
-Y t Nm 2 , where 



N 



4 8 

5 10 



(16) 



and m 2 = (m 2 _ i2 ,m 2 b ) T . The matrix N has eigenvector e.\ = (4,5) with eigenvalue 14, and 
another eigenvector with eigenvalue zero. 

We therefore consider what hierarchies are possible, starting with initial boundary con- 
ditions 



m H2 = \ m \ ■ ( 17 ) 



The suppression factor S is plotted in Fig. |3|. Several features are worthy of note. First, S 
is rather insensitive to 610, but is slightly improved for large particle content. In addition, 
even though the scalar masses are never pushed negative, S is degraded for very large 
h* = ht{M*). As above, this is caused by two- loop effects. Nevertheless, even for initial 
Yukawa couplings h* ~ 2, we see that \/S ~ 2 is possible. Scalars at 2 TeV again ameliorate 
a number of problems, and also stretch the discovery limits of planned future colliders. 

B. SU(5) 

Finally, we consider the minimal SU(5) model and extensions. The minimal SU(5) model 
has superpotential 

W = ^h t TTH u + V2h b TFH d + . . . , (18) 
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FIG. 3. The suppression factor S for SO(10) models evolving from M* = 2.4 x 10 18 GeV to 
Mq ~ 2 x 10 16 GeV with boundary conditions given by Eq. (|l7|). The parameter b\§ is the one-loop 
(3- function coefficient parametrizing this class of models, and h* is the value of the universal Yukawa 
coupling at the scale M*. 

where the quark and lepton fields are T (10) and F (5) and the Higgs fields are Ti u (5) and 

n d (5). 

We assume that h t and hb are equal and large, and that all other couplings are small 
relative to these. Since we take h t and to be near their quasi-fixed points, the latter 
restriction is not too severe.^] Note that we must take h t and to be large in order to 
suppress mn n and rriH d - This was an automatic and attractive feature of the minimal 
SO (10) models. 

The RGEs are given in the Appendix. In the minimal model, the /3-function coefficient 
is 65 = —3; it is larger in extensions of the minimal model, and so, as before, we take it to 
be a free parameter. 

The evolution matrix is 




6 Typically, the superpotential includes a term \TL u TfhLd, where S is the adjoint of SU(5). The 
requirement that colored Higgses be sufficiently heavy to prevent proton decay then implies A ~ g. 
However, in our scenarios, the first and second generation squarks are very massive, and so such 
constraints are relaxed by either a factor of (m^ght/m-heavy) 4 and an additional mixing angle or by 
a product of two additional mixing angles. 
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N 



3 6 

3 8 2 2 

4 4 4 

4 4 4 



(19) 



in the basis m 2 



m 



F ' 



m H d ) ' wnere Y 



Y t = Y\>. The leading eigenvector is 



e\ = (3,5,4,4), with eigenvalue 13. Suppression factors, with initial conditions 



3 2 
—rrin, 

4 rt « 



rrii 



m 



H d > 



(20) 



are plotted in Fig. ||. We find results similar to the SO (10) case. First and second gener- 
ation scalar masses as large as 2 TeV are allowed. As before, they significantly reduce the 
stringency of the experimental constraints on scalar degeneracy and CP- violating phases. 



V. CONCLUSIONS 

In this study we have examined the possibility that soft supersymmetry breaking scalar 
masses are of order mheavy at some high scale. The third generation scalar masses are 
driven to miight ~ TeV in the infrared by their large Yukawa couplings, while first and 
second generation scalar masses remain at mheavy This inverse hierarchy mechanism offers 
an appealing and natural way to satisfy the strong experimental constraints on first and 
second generation scalars in supersymmetric theories. 



13 



The third generation scalar masses are determined by infrared fixed points, so their 
precise values are insensitive to the initial conditions at, say, Mq (or to the exact value of 
Mq, for that matter). We have shown that with suitable assumptions and approximations, 
the scalar fixed points may be extracted from the RGEs analytically. We justified the validity 
of these approximations by more precise numerical calculations. 

For evolution below the GUT scale, we find that hierarchies of mheavy/^Hght ~ 4 may be 
created. The necessary boundary conditions are those of Eq. flTTp, with the Higgs squared 
masses double those of the squarks and sleptons. The simplicity of this boundary condition 
is in some sense a measure of the plausibility of the scenario. It is remarkable that such a 
simple boundary condition, perhaps the simplest imaginable next to complete scalar uni- 
versality, allows one to move the entire scalar mass scale up to 4 TeV without sacrificing 
naturalness. Note that for simplicity we have focussed on only the eigenvectors with largest 
eigenvalue. However, other eigenvectors with non-zero eigenvalue are also crunched, and so 
initial conditions in these directions may also be acceptable. It would be interesting to find 
theories that naturally generate such boundary conditions and to investigate the fixed point 
structure of other models. 

The scale mheavy ~ 4 TeV eliminates many dangerous supersymmetric effects; for exam- 
ple, supersymmetric contributions to EDMs, which in these scenarios scale like rrii ight /m^ ea:vy , 
may be easily below the experimental bounds even for 0(1) phases. However, these hier- 
archies are not sufficient to completely satisfy the most stringent experimental bounds by 
themselves. For example, for mheavy ~ 4 TeV, the K — K mass difference still requires 



rn 



i2/ m Lavy ~ 0.1, where m\ 2 is the off-diagonal entry of the squark mass matrix mixing the 
first and second generations. Nevertheless, such constraints are significantly weaker than in 
typical scenarios. 

For evolution above the GUT scale, the brief evolution interval makes it impossible to 
generate such large hierarchies. However, mheavy /flight ~ 2 is possible, and constraints on 
scalar degeneracy and CP- violating phases are still significantly weakened relative to the 
standard scenario in which all scalars are below 1 TeV. Note that this is opposite to the 
typical situation is unified models, where 0(1) generational splitting from unification-era 
evolution [IS] can lead to new contributions to flavor ||i~9| , p0|| and CP violation [21]. (A 



similar situation occurs in the presence of a right-handed neutrino with a large Yukawa 
coupling p2| .) 

Throughout this analysis, we have neglected off-diagonal scalar squared masses. For the 
squark masses, constraints on m 12 have already been noted above. The off-diagonal elements 
?7Ji3 an d rn 2 3 are also constrained: the requirement that the squark mass matrix be tachyon- 
free implies m^m^ < mi ight mheavy The off-diagonal masses are not driven to m lig h t and 
are, in fact, essentially RGE- invariant; the origins of their slight suppressions must therefore 
lie elsewhere. Depending on the extent of these suppressions, some interesting signals, for 
example, in the B system, may be observable in the near future [2j[[. 

Finally, we draw some implications for high energy colliders. In these scenarios, as in all 
scenarios with inverted scalar hierarchies, scalars with mass mheavy are beyond the reach of 
proposed colliders, such as the LHC and NLC, and can only be detected through their non- 
decoupling from super-oblique corrections |24[]. A goal of supersymmetric collider studies 
is to precisely measure superparticle properties, and thereby determine the supersymmetric 
parameters at the weak scale. It is then hoped that hints about physics at even smaller 
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distances can be found by evolving these parameters to very high energies. The fixed point 
analysis shows that in certain directions of parameter space, corresponding to eigenvectors 
of the scalar mass evolution matrix with large eigenvalues, even large masses at the high 
scale are exponentially suppressed when they are evolved to the weak scale. Accurate deter- 
mination of high scale boundary conditions in these directions therefore requires extremely 
precise measurements. 
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APPENDIX: 



In this appendix, we present the renormalization group equations for SU(5) and SO(10) 



GUTs above and below the GUT scale. As explained in Sec. |T|, gaugino masses and A 
terms are negligible in the scenarios we are considering. The relevant parameters are then 
the gauge couplings g, Yukawa couplings h, and scalar squared masses m 2 . The RGEs are 
conveniently expressed in terms of the variables 



t 



a 



Y 



In (M 2 /Q 2 ) 



g 1 



167T 2 

h 2 
16/T 2 



a 

47T 



X = m 2 x , 



(Al) 
(A2) 

(A3) 
(A4) 



where, as indicated in Eq. (|A4j), for notational simplicity we denote a scalar field's squared 
mass by the field itself. 

Schematically, the two-loop RGEs then have the form 



a = ba 
Y = Y 

X = X 



z2 



a 2 Y 



a 3 



-Y + a+Yl- 
-Y + Y 2 - aY 



aY 



a 2 



a 2 



(A5) 



where ( ) = 4. The constant b is the one-loop (5- function coefficient; the signs of all other 
terms are determined as indicated. Most of the two-loop terms are insignificant, as described 
in Sec. [IT]. Below, we present all one-loop terms and the dominant two-loop terms underlined 
above. 
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1. MSSM with Right-handed Neutrino 



The superpotential for the MSSM with a right-handed neutrino supermultiplet N is 

W = h t H u QU + h b H d QD + h T H d LE + h n H u LN . (A6) 
The RGEs for the three (GUT normalized) gauge couplings are 

& 3 = 3a 2 3 + a 2 3 {AY t + AY b ) 

<5 2 = -a\ + a 2 2 (QY t + QY b + 2Y T + 2Y n ) 

33. 2 _ 2 /26^ 14^ 18 6 \ 
«i = + a 2 y—Y t + + — Y T + -F n J . (A7) 

The Yukawa coupling RGEs are 

16 13 
Y t = Y t {-6Y t -Y b - Y n + -a 3 + 3a 2 + —a, 

+22Y t 2 + 5Y t Y b + 3Y t Y n + 5Y b 2 + y 6 y r + F T F n + 3Y 2 ] 

16 7 

n = Y b [-Y t - 6Y b -Y T + —a 3 + 3a 2 + —a x 

+5Y t 2 + 5F t n + Y t Y n + 22Y 2 + 3Y b Y T + 3Y 2 + Y T Y n ] 

Y T = Y T [-3Y b - AY T -Y n + 3a 2 + |di 

o 

+3Y t Y b + 3Y t Y n + 9Y 2 + 9Y b Y T + 10Y 2 + 3Y T Y n + 3Y 2 } 
Y n = Y n \-3Y t -Y T - AY n + 3a 2 + ^ 

+9Y 2 + 3Y t Y b + 9Y t Y n + 3y 6 Y T + 3Y 2 + 3Y T Y n + 10Y n 2 ] . (A8) 

Finally, assuming that the hypercharge trace tr [Y'm 2 ] vanishes, the scalar squared masses 
evolve as 

Q = -Y tSl - Y b s 2 + 10Y t 2 Sl + 10Y 6 2 s 2 + Y t Y n ( Sl + s 4 ) + Y b Y T {s 2 + s 3 ) 
f> = -2Y tSl + 1QY 2 S1 + 2Y t Y b ( Sl + s 2 ) + 2Y t Y n ( Sl + s 4 ) 

/j = -2y bS2 + i6F 6 2 s 2 + 2Y t Y b ( Sl + s 2 ) + 2y 6 y t ( S2 + s 3 ) 

L = -Y T s 3 - Y n s A + 6Y T 2 s 3 + QY 2 s 4 + 3Y t Y n ( Sl + s 4 ) + 3r 6 K(s 2 + s 3 ) 
£ = -2Y T s 3 + 8Y T 2 s 3 + 2Y r Y n (s 3 + s 4 ) + 6Y b Y T (s 2 + s 3 ) 
A> = -2Y n s 4 + 8Y 2 s 4 + 2Y T Y n (s 3 + s 4 ) + 6Y t Y n ( Sl + s 4 ) 
# u = -3y t si - Y n s 4 + 18y t 2 Sl + 3Y t Y b ( Sl + s 2 ) + 6Y; 2 s 4 + Y T Y n (s 3 + s 4 ) 

if d = -3y 6 s 2 - y r s 3 + i8>f s 2 + 3yy 6 ( Sl + S2 ) + 6y r 2 s 3 + y T y(s 3 + s 4 ) , (A9) 

where 

si = + Q + U 
s 2 = H d + Q + D 
s 3 = H d + L + E 

s A = H u + L + N . (A10) 

These RGEs are valid above the right-handed neutrino mass scale m^. The RGEs for 
the MSSM with minimal field content may be obtained, of course, by setting Y n = 0. 
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2. SO(10) above the GUT scale 



In the text, we consider a class of minimal SO(10) models with only one large Yukawa 
coupling. Here, for completeness, we present RGEs for the more general class of models 
specified by the superpotential 

W = htiHM2 + twWHi + --- , ( An ) 

where the quark and lepton fields are contained in if) (16), and H 2 (10) and H\ (10) are 
Higgs fields. In addition to the usual scalar mass terms, we also include the off-diagonal 
mass (H\H 2 + H 2 H\), which is allowed by all symmetries and is generated even if it 
is initially set to zero. The RGEs presented here correct certain one-loop terms given in 
Ref. ||19| , where the effects of this term were omitted. 

We consider models with arbitrary additional field content, but with interactions such 
that the dominant Yukawa couplings are h t and h b . The RGEs are then 

fi = -b w a 2 + a 2 (28Y t + 28Y b ) , (A12) 



Y t = Y t [-UY t - UY b + —a + 130Y, 2 + 2mY t Y b + 130Y 6 2 ] 

63 

Y b = Y b [-UY t - 14Y 6 +y & + 130Y t 2 + 26 °™ + 130Y 6 2 ] , (A13) 

and 

if; = -5Y tUl - 5Y b u 2 - 10Y tb H 12 

+90Y t 2 Ul + 90Y t Y b (ui + u 2 ) + 90Y b 2 u 2 + 180(Y t + Y b )Y tb H l2 
H 2 = -4Y tUl - AY tb H 12 

+mY 2 ui + 40Y t Y b ( Ul + u 2 ) + 120Y t Y tb H l2 + A0Y b Y tb H 12 
Hi = —4Y b u 2 - AY tb H X2 

+80Y 2 u 2 + 40Y t Y b (ui + u 2 ) + \2QY b Y tb Hi 2 + 40Y t Y tb H 12 
H 12 = -2Y tb {ui + u 2 ) - 2{Y t + Y b )H X2 

+20Y 4 F t6 (3 Ml + u 2 ) + 20Y 6 Y t6 (3M 2 + u x ) + (20Y, 2 + l2WY t Y b + 2VY b 2 )H X2 , (A14) 

where 

u x = 2iP + H 2 
u 2 = 2^j + H! 
H 12 = m 2 Hi2 

Y t l = y t y b , (M5) 

and 

b l0 = -1Q + S(R) . (A16) 

Common representations of SO(10) and their Dynkin indices are 10 (1), 16 (2), 45 (8), 54 
(12), 120 (28), and 126 (35). 
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In minimal SO(10) models, both up-type and down-type Higgs multiplets are contained 
in one SO(10) multiplet: TL 2 (10) D (H u , 7~Ld), where 7i u and H d are 5 and 5 representations 
of the SU(5) subgroup, respectively. The RGEs for these models are obtained by setting 
h b = Hi = TL\2 = in the equations above. The RGEs above, however, are also applicable 
to models in which the Higgs fields are contained in two separate SO(10) representations, 
with H 2 (10) D (H U ,H U ) and Hi (10) D (H d ,H d ). 

A useful consistency check follows from noting that, for h t = hb = hj\[2 and TCi = 7i 2 = 
Hi2 = TC/2, the two Higgs model reduces to a minimal model with Higgs field (Ti.2+'Hi) / 
and indeed, the resulting RGEs are identical to those for minimal models. 

3. SU(5) above the GUT scale 

We consider also the minimal SU(5) model and possible extensions. The minimal SU(5) 
model has superpotential 

W = -h t TTH u + V2h b TFH d + \iuH u H d + XHJM d + Ms trS 2 + \\' trS 3 , (A17) 
4 o 

where the quark and lepton fields are T (10) and F (5) and the Higgs fields are Ti u (5), H d 
(5), and S (24). 

The Yukawa couplings ht and hb are large and near their fixed point. We assume that all 
other couplings, such as A, A', fijj, and are small relative to these. The RGEs are then 

S = -b 5 a 2 + a 2 (l2Y t + UY b ) , (A18) 
96 

Y t = Y t [-9Y t - AY b + — a + 5AY t 2 + 2AY t Y b + 32Y b 2 ] 

Y b = Y h [-3Y t - 101^ + ^« + 18Y 2 + 30Y t Y b + QAY 2 } , (A19) 

5 

and 

t = -3Y t h - 2Y b t 2 + 3QY 2 h + 6Y t Y b (t! + t 2 ) + 32Y b % 
F = -AY b t 2 + \2Y t Y b {t x + t 2 ) + A8Y b % 
H u = -3Y t U + 36Y t % + UYtYbih + t 2 ) 

n d = -AY b t 2 + 12Y t Y b (ti + t 2 ) + 48Y b % , (A20) 

where 

h = 2T + H u 

t 2 = T + F + H d . (A21) 

In the minimal model, 65 = —3 in Eq. ( |A18[ ). For extensions of the minimal model, the 
above equations apply (again, assuming no large Yukawas), with the modification 

65 = -3 + S(R) , (A22) 

where S(R) is the sum of the Dynkin indices of the additional fields. For convenience, a 
few typical SU(5) representations and their Dynkin indices are listed here: 5 (|), 10 (§), 
15 (|), 24 (5), 45 (12), 50 (f ), and 75 (25). 
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